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Clustering of SZ clusters on a past light-cone: acoustic 
oscillations and constraints on dark energy 
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Abstract. We study the clustering of SZ-selected galaxy clusters on a past light-cone, particularly paying attention to the 
possibility of constraining properties of dark energy. The prospects of detecting baryonic features in the cluster power spectrum 
for a wide and shallow survey like PLANCK, and for an SPT-like narrow and deep survey are discussed. It is demonstrated 
that these future blank sky SZ surveys will have the capability to improve over the recently announced detection of baryonic 
oscillations based on the SDSS Luminous Red Galaxy (LRG) sample. We carry out parameter estimation using a Fisher matrix 
approach taking into account the anisotropic nature of the power spectrum due to redshift space and cosmological distortions. 
The clustering signal which is not too sensitive to systematic uncertainties serves as a valuable piece of information that in 
combination with other sources of data helps in breaking degeneracies between the cosmological parameters. 
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1. Introduction 

In the early 1970's it was recognized that acoustic waves in the 
radiation dominated matter prior to the epoch of recombina- 
tion of hydrogen in the Universe resulted in the characteristic 
pattern of maxima and minima in the post-recombination mat- 
ter power spectrum (Sunyaev & Zeldovich 1970; Peebles & Yu 
1970; Doroshkevich et al. 1978). These acoustic peaks depend 
on the size of the sound horizon and on the relative phases of 
the perturbations containing different masses at the moment of 
recombination. In the currently most favourable cosmological 
models with a Cold Dark Matter (CDM) component dominat- 
ing signifiantly over the baryonic part, the acoustic features in 
the matter power spectrum are strongly damped, reaching only 
~ 5% level for the "concordance" model (Bahcall et al. 1999; 
Spergel et al. 2003). For early description of the acoustic oscil- 
lations in the context of the CDM models see Blumenthal et al. 
(1988). For a given set of parameters of the Universe (Clbh 2 , 
Q„,/i 2 ) the position of the maxima and minima are fully deter- 
mined (e.g. Eisenstein & Hu 1998). Acoustic oscillations also 
leave their imprint on the angular perturbations of the Cosmic 
Microwave Background (CMB). MAXIMA- 1, 1 Boomerang, 2 
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WMAP, 3 VSA, 4 CBI, 5 and many other CMB experiments 
detected the first acoustic peaks in the CMB power spectrum 
with a high confidence level (Hanany et al. 2000; Netterfield 
et al. 2002; Bennett et al. 2003; Grainge et al. 2003; Pearson 
et al. 2003). These observations gave very important informa- 
tion about the key parameters of the Universe using the angular 
scale of acoustic features as rulers and taking into account the 
ratio of amplitudes of the different peaks. 

The existence of acoustic rulers in the Universe is of enor- 
mous importance since they permit us to measure the behavior 
of the Hubble parameter with redshift and also allow us to es- 
tablish the distance-redshift relation. This is especially impor- 
tant now when the discovery of dark energy (DE) is introducing 
more questions than the answers it provides. Different areas of 
observational cosmology (SNe la e.g. SCP , 6 high-z SN search; 
7 large scale structure surveys e.g. SDSS, 8 2dF 9 ; CMB exper- 
iments) provide evidence that the expansion of our Universe 
has been proceeding in an accelerated fashion since z ~ 0.75 
(Riess et al. 1998; Perlmutter et al. 1999; Spergel et al. 2003). 
Currently there is no physical understanding or even a reliable 
model for the DE. One of the first tasks to understand the na- 
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ture of DE will be the measurement of its equation of state 
parameter, w = P/p, and its possible evolution with time w(z). 
When we are equipped with a good standard ruler and are able 
to measure its angular behavior with redshift, we can obtain 
very valuable information about w(z). With the CMB data we 
can determine an angular diameter distance to the last scatter- 
ing surface with high precision. By combining this information 
with the measurement of acoustic peaks in the distribution of 
baryons at lower redshifts < z < 1 - 2, we will have unique 
information about the effective w and may even be able to deter- 
mine w(z) (Eisenstein & Hu 1998; Blake & Glazebrook 2003; 
Hu & Haiman 2003; Linder 2003; Seo & Eisenstein 2003). 

It was obvious since the first publications that acoustic os- 
cillations should also leave their imprint on the large scale 
structure of the Universe and thus influence e.g. the correla- 
tion function and the power spectrum of galaxies. The first suc- 
cessful detection of these features was presented by Eisenstein 
et al. (2005) who found traces of acoustic oscillations in the 
distribution of luminous red galaxies for which they had excel- 
lent measurements of angular positions and redshifts obtained 
by the SDSS collaboration. 

In this paper we discuss the opportunities that will be 
opened by the planned blank sky deep SZ cluster surveys which 
will be performed in the coming years. Clusters of galaxies 
are especially interesting objects for the study of acoustic fea- 
tures in the spatial distribution of objects since it has long been 
known that the clustering of clusters is an order of magnitude 
enhanced in comparison to galaxies (Bahcall & Soneira 1983; 
Kaiser 1984). Therefore even with smaller statistics of clusters 
it is possible to get useful results. 

In the next few years there will be very deep SZ cluster sur- 
veys of the restricted regions of the sky performed by several 
projects, e.g. APEX, 10 SZA, 11 AMI, 12 ACT 13 . Our analy- 
sis showed that the volume of these surveys and the number 
of possible cluster detections will be unfortunately insufficient 
for the search for the acoustic wiggles in the power spectrum. 
However, two planned surveys which will be carried out by the 
PLANCK Surveyor 14 spacecraft and the South Pole Telescope 
15 (SPT) have very good prospects for the detection of acous- 
tic features. PLANCK will make a shallow blank sky cluster 
survey permitting one to detect up to 20, 000 rich clusters of 
galaxies (e.g. Majumdar & Mohr 2004) with the bulk of ob- 
jects at z < 0.5, but will reach distances of z ~ 0.8. In contrast 
the SPT survey will observe deeper, but will cover only 10% 
of the sky. This survey is expected to detect up to 30, 000 clus- 
ters of galaxies (Majumdar & Mohr 2004), and many of them 
will be at significantly higher redshifts compared to the ones 
observed by PLANCK. Unfortunatelly it is not enough to mea- 
sure only the SZ flux or brightness of the clusters. In order to 
measure the equation of state of DE we need the redshift esti- 
mate for each cluster in the sample which will be hard and time 
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consuming work for many optical, X-ray and possibly radio 
astronomers. However, when this problem is solved, cosmol- 
ogists will have a unique sample of clusters of galaxies with 
good knowledge of their angular position, redshift, and hope- 
fully also mass. In this paper we investigate what limits to the 
DE equation of state might be obtained when these large exper- 
imental efforts are completed. It is obvious that in parallel other 
ways to measure w(z) will be implemented, but any additional 
and independent information will be useful. Especially impor- 
tant is that PLANCK and SPT surveys of clusters of galaxies 
will be performed in any case. Certainly, for many various pur- 
poses: (i) study of the redshift distribution of clusters, (ii) study 
of the properties of the clusters as a population, (iii) search 
for high-z clusters etc., we always need to estimate redshift. 
Therefore the information on the power spectrum of clusters 
of galaxies, acoustic wiggles and the subsequent determination 
of w is complementary but extremely important part of these 
surveys. 

The structure of this paper is as follows. In Sec. 2 we de- 
scribe an analytical model for a cluster power spectrum on a 
light-cone and calibrate it against the numerical simulations. 
Sec. 3 discusses the possibility of detecting baryonic oscilla- 
tions with the forthcoming SZ surveys. In Sec. 4 we carry out 
parameter forecasting using a Fisher matrix approach and Sec. 
5 contains our conclusions. 

2. Light-cone power spectrum of galaxy clusters 

In this section we present the theoretical model for calculat- 
ing the cluster power spectrum on our past light-cone. In order 
to assess the accuracy of the theoretical description, we make 
comparisons with the VIRGO Hubble Volume simulation out- 
puts. We start with a very brief description of the VIRGO simu- 
lations and proceed with the calculation of power spectra using 
light-cone cluster catalogs provided by the VIRGO Consortium 

16 

2.1. Cluster power spectra from VIRGO simulations 

We use outputs from the ACDM Hubble Volume simulation 
that was run in a 3000 hr x Mpc comoving box with a parti- 
cle mass of 2.25 • 1Q 12 h~ l M Q . The other simulation parameters 
were as follows: Q m = 0.3, 0\ = 0.7, Q. b = 0.04, h = 0.7 and 
erg = 0.9 (for further details see e.g. Evrard et al. 2002). In our 
study we used z = snapshot, SphereB and OctantB light-cone 
catalogs 17 . The z = cluster catalog used a friend-of-friend 
scheme with a linking length b = 0.164 for cluster identifica- 
tion, while for the light-cone outputs the spherical overdensity 
method with the overdensity 200 relative to the critical density 
was applied. The minimum number of particles per cluster is 
30 and 12 for the z = and light-cone catalogs, respectively. 

To calculate the power spectrum we follow the direct 
method of Feldman et al. (1994) (FKP), which is shown to be 
optimal for sufficiently large A:-modes i.e. k » 1/L, where L 
is the typical spatial extent of a survey volume (Tegmark et al. 

16 http ://w w w. mpa-garching . mpg.de/Virgo/ 

17 For the exact description of these catalogs see Evrard et al. (2002) 
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1998). Because FFTs are used to achieve significant speedup 
for Fourier sum calculations, we first have to find the density 
field on a grid. To this end we use the Triangular Shaped Cloud 
(TSC) (Hockney & Eastwood 1988) mass assignment scheme. 
Thus our density field is a filtered version of the underlying 
field, and as shown in Jing (2005), the real power spectrum P 
can be expressed as the following sum over aliases (correct 
again for the case k » 1/L): 



Praw(k) * J] S 2(k ' } Z + 2fe Nn)P(k + 2fc N n) + 

i^nV 2 (k + 2fe N n), (1) 



neZ 3 



neZ 3 



where the raw power spectrum: 

Praw(k) = <|5 g (k)| 2 > 



and Fourier transform of the overdensity field on a grid is cal- 
culated as usual: 



W = ^2[«g(rg))-«S(r g )le^ k . 



(3) 



Here n g is the number density field on a grid without any selec- 
tion effect corrections, S is a selection function that also incor- 
porates survey geometry (i.e. 5=0 outside of survey bound- 
aries), h is the mean underlying number density and the sum 
runs over all grid cells. A Fourier transform of the selection 
function S (k) in Eq. (1) is normalized so that 5(0) = 1 and the 
mass assignment window in the case of the TSC scheme can be 
expressed as: 



nV(k) 



Ilsin(f|)l 3 



n ^ 

1.1 2* N 



(4) 



The second term on the left-hand side of Eq.(l) is the shot noise 
contribution and in the case of the TSC filter can be shown to 
give the following result (Jing 2005): 



SN 



1 - sin 



21c N 



2 , nk 
+ — sin^ — — 
15 \2k N 



(5) 



Summarized very briefly, our power spectrum calculation con- 
sists of the following steps: 

1 . Determination of the selection function S (including survey 
geometry) and mean underlying number density n, 

2. Calculation of the overdensity field on a grid using the TSC 
mass assignment scheme and its Fourier transform as given 
in Eq. (3), 

3. Subtraction of the shot noise term (Eq. (5)) from the raw 
power spectrum (Eq. (2)), 

4. Isotropization of the shot noise corrected power spectrum, 
i.e. averaging over fc-space shells, 

5. Application of normalization correction due to selection ef- 
fects, i.e. dividing by £ k 5 2 (k), 

6. Deconvolving the smearing effect of the TSC mass assign- 
ment. 




(2) Fig- 1- Mass conversion from M^ m to M 180 for various redshifts. 



The "sharpening" in the last step is done using an iterative 
method as described in Jing (2005) with the only difference that 
here we do not approximate the power spectrum simply with a 
power law, but also allow for a running of the spectral index, 
i.e. we approximate it with a parabola in log-log coordinates. 

When calculating the power spectrum in the above de- 
scribed way we assumed that the influence of selection/survey 
geometry effects on the power spectrum are separable. This is 
not the case for large scales (k ~ 1 /L) and also if too narrow 
(i.e. Ak < 1/L) power spectrum bins are used. In the following 
we always make the power spectrum binning broad enough so 
that the neighboring bins can be safely assumed to be uncorre- 
cted. 

The power spectrum error is estimated using the simple 
"mode counting" result of FKP (see also Tegmark et al. 1998): 



AP _ f 



VeffVk' 



(6) 



where Vk = 4-7Tk 2 Ak/(2n) 3 is the volume of the k- space shell 
and V ei f is the effective volume given by: 



eff 



(7) 



Here dV c is a comoving volume element and the weight func- 
tion: 



W(z)< 



const 

n(z) 

l+«(z)P 



for volume weighting 
for number weighting 
for an optimal FKP weighting. 



(8) 



The halos of the VIRGO Hubble Volume simulations were 
identified using the spherical overdensity algorithm with an 
overdensity of 200 with respect to the critical density at the 
identification epoch (Evrard et al. 2002) (we denote the cor- 
responding mass Mjqq), while detailed comparisons suggest 
that the Press-Schechter type of analytical calculations (Press 
& Schechter 1974; Mo & White 1996; Shefh & Tormen 1999; 
Sheth et al. 2001) provide a good match to simulations if an 
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overdensity 180 with respect to the background density is used 
(corresponding mass Migo) (e.g. Jenkins et al. 2001). In order 
to convert from one mass definition to the other we assume that 
the density profile of clusters is given by the NFW (Navarro 
et al. 1997) profile: 



P(0 



(9) 



and the concentration parameter c = r v /r s (r v -virial radius) and 
its evolution as a function of virial mass M v is given as follows 
(Bullock et al. 2001): 



c(M v ) 



M v 



l+z\Af„(z = 0) 



-0.13 



(10) 



Here M, is a standard nonlinear mass scale defined through 
cr(M t ,z) = S c (z), where cr 2 (M,z) is the variance of the lin- 
early evolved density field on the comoving scale correspond- 
ing to the mass M at redshift z and d c (z) is the spherical collapse 
threshold e.g. in Einstein-de Sitter model 6 C (0) = 1.686. 
Then the mass within the radius r can be expressed as: 



M(< r) = 4 W 3 /(-^) = — AQ mPc , 
where 



fix) 



In 1 + 



1 



1 



1 + x 



(11) 



(12) 



p c is the critical density and A is the halo overdensity with re- 
spect to the background matter density at the epoch of halo 
identification. 

In order to convert halo mass M corresponding to the over- 
density A to the one corresponding to the overdensity A' we 
proceed as follows: 

1. From Eq. (11) determine the radius r corresponding to the 
mass M and overdensity A, 

2. Solve 



Av/ 



c(r v )r 



-A/ 



1 



c(r v ) 



= 



(13) 



for virial radius r v . Here A v is the virial overdensity that 
we find numerically solving the spherical tophat collapse 
model (fitting formulae for A v for some cosmological mod- 
els are given in Bryan & Norman (1998)), 
3. Solve 



Av/f-^WyU-U 

\c(r v y) \c(r v )j 



(14) 



for r', 

4. From Eq. (11) find M' corresponding to the radius r' and 
overdensity A' . 

The results of this mass conversion from M^ 00 to M^o 
are shown in Fig. 1 as a fractional increase in mass / = 
A^i8o/A^20o ~~ 1 f° r different redshifts. We see that especially for 
low redshift clusters this mass change can reach up to 50%. In 
Fig. 2 we demonstrate the importance of the mass conversion in 
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0.0 < z < 0.2 (x 1000) 
0.0 < z < 0.5 (x 100) 
0.0<z<1.0 (x 10) 
0.0<z< 1.45 



+++ + ■> 



M 1M (h-MJ 



Fig. 2. Light-cone mass functions for different redshift intervals. For 
clarity the curves have been shifted by the factors given in the legend. 
Boxes show simulation results with the applied mass conversion from 
M^qq to Mi go whereas crosses are the results without any conversion. 



order to get agreement with the analytical mass function calcu- 
lations. Here the mass function on a light-cone was calculated 
as follows: 



n LC (M m , < z) = - 



Jn(M m ,z)^dz 



(15) 



Here dV c is a comoving volume element and n(M,z) is a 
mass function as described in Sheth & Tormen (1999), which 
is known to give a very good description of N-body results 
(Jenkins et al. 2001). 

In the upper panel of Fig. 3 we show various correc- 
tions needed to achieve a reliable estimate of the power spec- 
trum of the underlying cluster distribution while the lower 
panel demonstrates the consistency of the applied "sharpen- 
ing" scheme. Here we have used the z = cluster catalog to 
allow for a comparison with the results presented in Colberg 
et al. (2000). The lower mass for the cluster selection was 
taken to be 1.0 • 10 14 /z~'M o in order to get the total number 
of objects equal to ~ 915,000 as was used in Colberg et al. 
(2000). Also we have selected a spherical volume out of the 
full box to test how well the geometry correction works. The 
results of this comparison are given in Fig. 4. We see that the 
Colberg et al. (2000) power spectrum agrees with our calcu- 
lations at the largest scales; however, for the smaller scales it 
drops below our results. We suspect that their correction for the 
grid smoothing effect was insufficient, although in their paper 
they do not describe how the power spectrum was calculated. 
As can be seen from the figure the shape of our cluster power 
spectrum agrees very well with the linear theory matter power 
spectrum up to the scale k ~ 0. 15 h Mpc~' . Clearly with such a 
large number of clusters (~ 477, 000 inside our spherical vol- 
ume) baryonic oscillations are easily detectable and the cor- 
responding "smooth" model without them is disfavored. The 
theoretical matter power spectra were calculated as described 



G. Hutsi: Clustering of SZ clusters 



5 




M 18() > 10 14 h'' M Q , 128 3 TSC-grid, z = 



k |h Mpc' 



64 TSC-grid, corrected 

61 TSC-grid 

128^ TSC-grid. corrected 

128', TSC-grid 

256, TSC-grid, corrected 

256 TSC-grid 




Fig. 3. Upper panel: various corrections applied to reach the final 
power spectrum estimate. Lower panel: results of the consistency test 
for the "sharpening" scheme using different grid sizes. 



in Eisenstein & Hu (1998). Using the z = cluster catalog 
we also calculate power spectra and two-point correlation func- 
tions for various lower mass cutoffs. These results are presented 
in Fig. 5 where the left-hand panels show power spectra divided 
by the smooth model without baryonic oscillations and right- 
hand panels the respective correlation functions. Here the up- 
permost power spectrum is the same as the one given in Fig. 4. 
Solid/dotted lines show theoretical models with/without bary- 
onic oscillations. Correlation functions were calculated using 
the estimator given by Landy & Szalay (1993): 



DD - 2DR + RR 
RR ' 



(16) 



which has minimal variance for a Poisson process. Here DD, 
DR and RR represent the respective normalized data-data, data- 
random and random-random pair counts in a given distance 
range. Random catalogs were generated with ten times the 
number of objects in the main catalogs. The survey geome- 
try was again taken to be a spherical volume reaching redshift 
z — 0.58. The number of objects corresponding to the lower 
mass cutoffs of 1.0 • 10 14 , 2.0 • 10 14 , 3.0 • 10 14 , 4.0 • 10 14 and 
5.0 • 10 14 hr l M B (friend-of-friend masses) were respectively 



256 J TSC-grid 1 
theoretical ps 
theoretical 'smooth' ps 
Colberg et al. 2000 




k [h Mpc' 1 ] 

Fig. 4. Upper panel: Power spectrum of clusters more massive than 
1.0 • 10 14 /r'M o from z = simulation box. Points with errorbars 
show the power spectrum of clusters inside a spherical subvolume 
with comoving radius of 1500 h' 1 Mpc extracted from the full sim- 
ulation box. Stars present results obtained by Colberg et al. (2000) 
and solid/dashed lines are model power spectra with/without acoustic 
oscillations. Lower panel: as above, except all the curves have been 
divided by the "smooth" model without acoustic oscillations. 



476, 634, 167, 898, 80, 81 1, 45, 836 and 27, 955. For the corre- 
lation function we have shown only a simple Poissonian errors: 



1+f 



(17) 



(Errors due to the terms DR and RR can be neglected because 
of the much larger number of available pairs.) These errorbars 
are an underestimate of the true variance, being only exact for 
the Poissonian point process. To estimate true errors, one needs 
to have a knowledge about the 3- and 4-point correlation func- 
tion of the clustering pattern. We are not going to elaborate 
further on these issues since the following analysis is based 
solely on the power spectrum. From Fig. 5 we see that bary- 
onic features are visible down to the case with the lowest num- 
ber of clusters. The only exception is the correlation function 
with the lowest number of objects which shows rather noisy 
behavior near the expected baryonic bump. For the least mas- 
sive systems, due to the nonlinear evolution, the power spec- 
trum starts to rise at k > 0.15 /zMpc -1 . For the more massive 
systems, on the other hand, the opposite trend is visible i.e. a 
decrease of power. This is due to the cluster formation which 
can be viewed as a smoothing filter acting on an initial density 
perturbation field. 
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0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2 60 80 100 120 140 

k|hMpc"'| rlh^Mpc] 



Fig. 5. Power spectra (left panels) and correlation functions (right panels) for the lower mass cutoffs of 1.0 • 10 14 , 2.0 • 10 14 , 3.0 • 10 14 , 4.0 • 10 14 
and 5.0 • 10 14 h~ l M e (friend-of-friend masses). Power spectra have been divided by the model spectra without acoustic features. The number 
of objects inside a spherical survey volume reaching z = 0.58 is given in each panel. Solid/dotted lines show theoretical models with/without 
baryonic oscillations. 



2.2. Comparison with the analytical description: 
accuracy of the biasing scheme 

The analytical power spectrum of clusters with masses M > 
Mi ow on a light-cone P^ c (k; > Mi ow ) is calculated as presented 
in Yamamoto et al. (1999) with a slight modification to allow 
for various weight functions, so: 



P\ c (k;> M low ) = 



Mnax 

J W 2 (z)P c (k;>M low ,z)^dz 

Zmax 

S ^(z)f dz 



(18) 



where the weight function W(z) is given in Eq.(8). There the 
number density of objects is provided by the cumulative mass 
function at redshift z: 



n(> M low ,z) 



j n(M, 



z)dM. 



(19) 



The power spectrum of clusters more massive than M\ ovl at red- 
shift z is given as: 

P c (k; > M low , Z ) = D\ (z)b 2 eff (> M low ,z)P(k, z = 0), (20) 
where the effective bias parameter: 



b e s(> Mi ow ,z) = 



f b(M,z)n(M,z)dM 
n(>M low ,z) 



(21) 
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D+(z) is the growing mode of linear density fluctuations nor- 
malized such that D + {z = 0) = 1 and P(k,z = 0) is the matter 
power spectrum at the current epoch, which is calculated using 
the transfer functions presented in Eisenstein & Hu (1998). 

For the mass function n(M, z) and bias parameter b(M, z) we 
use both Press-Schechter (PS) (Press & Schechter 1974) and 
Sheth-Tormen (ST) (Sheth & Tormen 1999) prescriptions. It is 
well known that PS mass function underpredicts the number 
density of massive objects (Sheth & Tormen 1999). PS overes- 
timates while ST underestimates the bias parameter for massive 
halos (especially at larger redshifts) (Sheth & Tormen 1999). 
PS underestimation of number density turns out to be approx- 
imately compensated for by its overestimation of bias param- 
eter, and as such, we get the best agreement with the numer- 
ical light-cone power spectra using a plain PS approach. This 
is demonstrated in Fig. 6 where we show the light-cone power 
spectra for various values of lower mass cutoff M low . Results in 
the upper panel apply for the SphereB cluster catalog (reaching 
redshift z ~ 0.58) whereas the ones on the lower part of the 
figure correspond to the OctantB catalog (reaching z ~ 1 -46) of 
VIRGO Hubble Volume simulation outputs. We obtained the 
best agreement if clusters were selected using M\%q, but in bias 
calculations virial mass M v ; r was used instead. 

Overall the agreement between the numerical results and 
an analytical description is better than 20%. 



3. SZ clusters and baryonic oscillations 

The "concordance" cosmological model predicts oscillations in 
the matter power spectrum with a relative amplitude of ~ 5% 
(see Fig. 7). The correspondance of the peaks in the matter 
power spectrum to the ones in the CMB angular power spec- 
trum is also given in Fig. 7 (The first vertical line represents 
the position of the 1st CMB acoustic peak etc.). We see that 
on small scales the corresponding features are out of phase. 
This is due to the so-called velocity overshoot, meaning that at 
those scales the growing mode of density fluctuations is mostly 
sourced by the velocity perturbations (Sunyaev & Zeldovich 
1970; Eisenstein & Hu 1998). At larger scales, on the other 
hand, fluctuations in the density provide a dominant source 
term and so the corresponding features in the power spectra are 
in phase. Each rise and fall in the matter power spectrum corre- 
sponds to twice as many features in the CMB angular spectrum. 
This is a generic property of the models with a dominating 
CDM component. Purely baryonic models on the other hand 
would have oscillations with the same frequency as in the an- 
gular spectrum of the CMB. 18 In order to see these features in 
the matter power spectrum one needs a "tracer" population of 
objects whose clustering properties with respect to the under- 
lying dark matter distribution are reasonably well understood. 
These objects should have high enough number density to re- 
duce discreteness noise on one hand, but on the other hand they 
should fill as large of a comoving volume as possible to de- 
crease cosmic variance. In general one wants to find an optimal 



s 




Fig. 6. Upper panel: light-cone power spectra for SphereB output of 
VIRGO ACDM Hubble Volume simulation for different lower mass 
cutoffs (see legend). Lower panel: same as above but for OctantB out- 
put. 



solution of these two degrees of freedom 19 , as to maximize the 
obtainable effective volume (see Eq. (7)) for a fixed observa- 
tional effort. Currently the largest effective volume amongst all 
of the available surveys is provided by the SDSS LRG sample 
(Eisenstein et al. 2005). The analysis of this sample yielded a 
detection of a clear acoustic feature in the spatial two point 
correlation function. Future projects such as the K.A.O.S. 20 
galaxy redshift survey has as one of its main scientific targets 
the detection of baryon oscillations in the spatial clustering of 
high-z galaxies. The possibility of using the aforementioned 
galaxy redshift surveys to measure the sound horizon has been 
discussed in several papers e.g. Blake & Glazebrook (2003); 
Linder (2003); Seo & Eisenstein (2003); Hu & Haiman (2003) 
and a similar discussion in the context of photometric redshift 
surveys is given in Blake & Bridle (2005) (see also the discus- 
sion in Seo & Eisenstein 2003). 



http ://cmb. as . arizona. edu/~ ei sens te/acousticpeak/ 



19 One also has to consider the steepness of the luminosity function 
of those objects. 

20 http://www.noao.edu/kaos/ 
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Fig. 7. Theoretical dark matter power spectra with a smooth compo- 
nent divided out (using transfer functions as given by Eisenstein & 
Hu (1998)) for the WMAP "concordance" model (solid line) and for 
the model used in VIRGO Hubble Volume simulations (dashed line). 
The numbered vertical lines show the locations of the corresponding 
peaks in the CMB angular power spectrum e.g. 1. corresponds to the 
1st acoustic peak etc. 

Here instead of galaxies we discuss the possibility of using 
SZ-selected galaxy clusters for that purpose. Some calculations 
related to the SPT-type of SZ survey were also presented in Hu 
& Haiman (2003). It is clear from Figs. 4 and 7 that with wide 
field galaxy cluster surveys we should be especially sensitive 
to the scales that correspond to the 2nd and 3rd acoustic peaks 
in the CMB angular power spectrum. 

A few advantages of using galaxy clusters compared to the 
galaxies are: 

- With relatively small cluster samples it is possible to probe 
large cosmological volumes (thus reducing cosmic vari- 
ance). 

- The clustering signal of galaxy clusters is amplified with 
respect that of galaxies. 

- The relation with respect to the underlying dark matter field 
is rather well understood and also redshift space distortions 
are managable since "fingers of god" could be avoided. 21 

The biggest disadvantage is a rather low number density i.e. 
high shot noise contribution. 

3.1. SZ-selected clusters. Mass-observable relations 

In order to compare observations with the models one has to 
establish mass-observable relations and also specify survey se- 

21 This is only true when we have enough galaxy redshifts per clus- 
ter, so that one can average down to something close to the center of 
mass velocity. One might also try to exploit the fact that the bright 
central cD galaxies have small velocities with respect to the rest of 
the cluster (see the simulation results by Berlind et al. (2003)). Even 
if one is able to find a good estimate for the center of mass velocity 
for each cluster in the sample, the redshift space distortions on reason- 
ably large scales would still deviate from the simple linear prediction 
(Scoccimarro 2004). 



lection criteria. Here for the sake of simplicity we assume that 
all the clusters remain unresolved i.e. we assume that our sam- 
ple is effectively flux-selected. This is a rather good approxi- 
mation for the case of PLANCK, but for surveys like ACT and 
SPT extra complications will arise since part of the cluster pop- 
ulation will be resolved and so the selection function has one 
additional degree of freedom, namely surface brightness. As 
our aim here is not to give any detailed predictions for a par- 
ticular survey these assumptions seem to be quite reasonable. 
The change in detected flux towards a galaxy cluster due to the 
thermal SZ effect can be expressed as: 



F(M,z,x) 



ho-T g(x)f b M-kT(M,z) 



jj. e m v m e c l 



d A (z) 2 



(22) 



where / = 2(kT cmb ) 3 /(he) 2 2.7 • 10 n mJy/sr, f h is the cluster 
baryonic fraction which we take to be equal to the cosmic av- 
erage Ob/^m, ^a(z) is angular diameter distance to the cluster, 
/j. e = 2/(1 + X) for the case of fully ionized plasma with neg- 
ligible metallicity (we take /i e = 1.14). The spectral function 
g(x) is given as follows (Sunyaev & Zeldovich 1980): 



*(*) = 



xcoth(-) -4 
2 



(e x - l) 2 ' 



(23) 



where the dimensionless frequency x = hv/kT cm b 0.0176 • 
v(GHz). For the mass-temperature relation we assume a simple 
virial scaling (Bryan & Norman 1998): 



kT(M,z)=A-[A c (z)E(z) 2 ]~Ml 



(24) 



Here A c (z) is a critical collapse overdensity with respect to the 
critical density at redshift z and E(z) = H{z)/Hq. The normal- 
izing constant A is determined so as to obtain a good match for 
the SZ cluster number counts from the state-of-the-art hydrody- 
namical simulations by White et al. (2002). These simulations 
included gas cooling processes and also feedback from super- 
novae and galactic winds. If we measure kT in keV and M in 
units of /i _1 M , then a good fit can be obtained if A 1.0- 10~ 10 
as seen in Fig. 8. Here we present results both for ST and PS 
mass functions. 

In reality the mass-observable relations are currently rather 
poorly known but one may argue that planned surveys with the 
yields of tens of thousands of galaxy clusters have significant 
power for "self-calibration" (Majumdar & Mohr 2004). As also 
shown in Majumdar & Mohr (2004) a much better approach 
would be to establish these scaling laws using external mass 
determinations (e.g. through lensing studies) for a subset of a 
complete sample. For a more precise modeling of the selection 
effects one also has to consider scatter around these mean rela- 
tions. These issues can be settled once we obtain a real sample. 
Moreover, the clustering as compared to the number count of 
objects is much less sensitive to the uncertainties in the precise 
knowledge of the selection effects. Here the selection effects 
enter when relating the clustering of tracer objects to the un- 
derlying dark matter i.e. while determining the effective bias of 
objects. As it turns out (see Sec. 4), future large cluster sam- 
ples are able to provide a good estimate of the effective bias 
themselves through the redshift space distortions. 
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Fig. 8. Cumulative number counts of SZ clusters for observing fre- 
quency 150 GHz. Crosses show simulation results by White et al. 
(2002) whereas solid/dashed lines correspond to analytical results as- 
suming a ST/PS mass function. 

3.2. Accuracy of the power spectrum determination 

We investigate how well we can determine the power spec- 
trum with SZ surveys having various sensitivity limits, specif- 
ically concentrating on the range 1 ... 25 mJy. Our results for 
AP/P (or equivalently for the effective volume, see Eq. (6)) 
and for the number of detectable clusters are given in Fig. 9. 
Here the upper panel assumes full sky coverage and the lower 
one applies for one octant of the sky. With the solid lines we 
have plotted the fractional accuracy AP/ P achieved for differ- 
ent lower flux limits and follow-up survey depths applying the 
FKP weighting scheme. The wavenumber k in the calculations 
was taken to be 0.05 /zMpc -1 , which is close to the first ma- 
jor acoustic feature in the expected matter power spectrum (see 
Fig. 7). The bin width Ak = 0.005 h Mpc~' is large enough so 
that even for the shallowest surveys reaching only z ~ 0.2 the 
power spectrum bins can be assumed to be independent. Here 
and also in the following we perform analytical calculations 
for the observational frequency 150 GHz and assume cosmo- 
logical parameters consistent with the WMAP "concordance" 
model (Spergel et al. 2003). 

Light-cone power spectra for SZ clusters are calculated us- 
ing Eq. (18) where the lower integration boundary M\ ovl {F, z) 
is given by Eq. (22) and (24). Then AP/P is found using Eq. 
(6), (7) and (8) and the mean underlying number density h(z) is 
given by Eq. (19). Also we have taken into account the increase 
of the isotropized power spectrum due to linear redshift space 
distortions by a factor of 1 + 2/3/3 + /3 2 /5 (Kaiser 1987), where 
r - _L . dlnD + 

" b cS d In a ' 

The solid lines in Fig. 9 starting from below correspond to 
AP/P values of 6%, 5%, 4%, ... in the upper panel and 10%, 
9%, 8%, ... in the lower one. With the dashed lines we have 
plotted the number of clusters. Moving from the lower right 
to the upper left each line represents a factor of two increase 
in number with the starting values being 20,000 and 10,000 
in upper and lower panels, respectively. The flattening out of 
AP/P and cluster number curves at low fluxes is due to the 
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Fig. 9. Upper panel: AP/P (given in %) and number of clusters for 
flux-limited surveys with various sensitivities and survey depths at 150 
GHz assuming full sky coverage. Solid lines show AP/P isocontours 
whereas constant number contours are given with dashed lines. Lower 
panel: the same as above only for one octant of the sky. 

imposed lower mass cutoff 1.0 • \0 u hT l M Q . Thus below some 
flux limit we see all the clusters inside a specified volume that 
have masses above that cutoff value. 

In general, by increasing the volume of the survey we also 
boost the shot noise contribution due to the decreasing num- 
ber density of distant objects. As seen from the figure- in the 
case of FKP weighting- the accuracy of power spectrum esti- 
mate does not degrade as we move to further distances since we 
down weight the contribution of the far away objects in such a 
way as to compensate for the increase in shot noise. For too 
small survey volumes, on the other hand, limits on achievable 
accuracy are set by the growing importance of the cosmic vari- 
ance. The FKP weighting scheme is not strictly optimal in our 
case since it was derived assuming a fixed i.e. non-evolving 
underlying power spectrum. Certainly, for a better scheme one 
should weight down the contribution of the far away objects 
slightly more mildly since the clustering strength of these ob- 
jects is higher. This type of weighting method, which is able 
to handle at least the case with an evolving amplitude, is pre- 
sented in Percival et al. (2004). Nevertheless in the following 
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calculations for simplicity we still apply the FKP weight func- 
tion. 

The results in Fig. 9 assumed that we have a full follow-up 
such that we are able to obtain all the redshifts of the clus- 
ters detected by an imaging survey. Also the sample was as- 
sumed to be purely flux-selected which is rather unrealistic for 
real experiments e.g. for PLANCK many clusters remain unde- 
tected due to rather poor angular resolution or oppositely in the 
case of SPT some fraction of clusters will be "resolved out" 
and a significant amount of signal will be lost. Many objects 
might remain undetected for these reasons. If the systems that 
are left out are low mass clusters (as in the case of PLANCK) 
then our power spectrum estimate might actually be almost as 
good as before since low mass systems are relatively weakly 
clustered and as such they do not contribute significantly to 
the total signal. This can be seen in Fig. 10 where in the up- 
per panel we have shown the influence of changing the lower 
mass cutoff Mi ow on AP/P. On the lower panel the respective 
number of clusters is given. These calculations are done for 
two different survey types: (1) solid lines represent results for 
a shallow SZ survey covering the full sky and reaching red- 
shift Z]im = 0.6 with a sensitivity limit Fi ow = 17 mJy at 150 
GHz, (2) dashed lines represent a deep survey (with no upper 
redshift cutoff imposed) covering 1/8 of the sky with a flux 
limit Fi ow = 5 mJy at 150 GHz. The first might be applica- 
ble for the case of the PLANCK mission 22 and the second for 
the SPT cluster survey. In practice the measurements will be 
performed in many frequency channels which helps to sepa- 
rate clusters from other foreground sources due to the specific 
frequency behavior of the thermal SZ effect. Here for simplic- 
ity we have chosen the sensitivity limits corresponding to the 
"weakest" of the channels available for SZ purposes. We see 
that for a full sky with the 20, 000 most massive clusters up 
to redshift z ~ 0.6, one could obtain an estimate of the power 
spectrum at k = 0.05 h Mpc~' with a fractional error below 5% 
while for the one octant of the sky with less than ~ 25, 000 
clusters we always stay above 5% accuracy. 

The previously described approach where we use all the 
data to obtain a single combined estimate of the power spec- 
trum is appropriate if we only intend to place constraints on 
n DE . Combining this power spectrum estimate which is sen- 
sitive to Q. m h with the CMB constraint on Q. m h 2 gives us Q. m 
and h separately. Additionally, knowing the geometry of the 
Universe from CMB measurements gives us immediately an 
estimate of £2de- Since in the majority of the DE models DE 
starts to dominate relatively recently, driving the Universe fur- 
thest from the plain Einstein-de Sitter behavior, the best red- 
shift to complement the CMB data is at z — 0. 

On the other hand if our aim is to constrain the equation 
of state parameter w and its possible change in time it is es- 
sential to measure the power spectrum at different redshifts. 
This leads to the question of how to bin up the sample in 
redshift? Certainly there are optimal ways of combining data, 
but unfortunately they all depend on the way we choose to 



22 Using the spectral dependence of thermal SZ effect (see Eq. (23)) 
we see that the 30 mJy sensitivity of PLANCK's 353 GHz channel 
corresponds to 17 mJy at 153 GHz. 
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Fig. 10. AP/P (upper panel) and number of clusters (lower panel) as 
a function of the lower mass cutoff M low . Solid lines correspond to 
the full sky shallow SZ survey and dashed lines to the deep survey 
covering one octant of the sky. 

parametrize our model for DE. Recently Huterer & Starkman 
(2003) argued that in the absence of a theoretically well moti- 
vated parametrization one should use a stepwise function with 
the value W{ in the z'-th redshift bin and let the data itself deter- 
mine which combinations of w>\ will be well constrained. One 
can then reconstruct the behavior of w with the redshift as a lin- 
ear combination of the "cleanest" eigenmodes. Here we are not 
trying to implement that kind of general parametrization since 
as a first step it should be sufficient to determine an effective 
constant w and see whether it deviates from the currently most 
well motivated w = - 1 . Therefore, in the following we mostly 
investigate the case with a constant equation of state parameter 
wo. The redshift binning is chosen so as to get equal relative ac- 
curacies of the power spectrum in each bin. The results of this 
binning procedure for the above described two types of survey 
are given in Fig. 1 1 with the solid lines corresponding to the 
shallow one. For the shallow (deep) survey we have assumed 3 
(4) redshift bins. In the inset the upper curves show the relative 
accuracy achievable in each redshift bin while the lower lines 
correspond to the full sample without any binning and so co- 
incide with the lines shown in the upper panel of Fig. 10. This 
kind of redshift division is also used in the following parameter 
estimation section. 
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Fig. 11. "Equal accuracy" redshift binning for different lower mass 
cutoffs. Solid (dashed) lines correspond to the shallow (deep) SZ sur- 
vey with 3 (4) redshift bins. In the inset the upper curves show the 
relative accuracy achievable in each redshift bin while the lower lines 
correspond to the full sample without any binning. 



3.3. Prospects of detecting baryonic "wiggles". 
Comparison with SDSS LRG 

The relative amplitude of acoustic oscillations in the mater 
power spectrum for the WMAP "concordance" model as well 
as in the power spectrum used in VIRGO Hubble Volume sim- 
ulations is shown in Fig. 7. Here the smooth component was 
divided out using the fitting formulae provided by Eisenstein 
& Hu (1998). In order to be able to detect these features the ac- 
curacy of the power spectrum determination should be of com- 
parable size i.e. AP/P ~ 5%. This implies that for a full sky 
survey one needs on the order of 25,000 galaxy clusters in- 
side the volume with limiting redshift z ~ 0.6 as seen from the 
upper panel of Fig. 10. It is clear that with clusters one might 
hope to detect only the first few acoustic signatures, e.g. the 
major features at k ~ 0.045 and ~ 0.075 h Mpc -1 , since they 
are too rare objects to enable the sampling of the smaller scale 
density field. 

Since the study of the SDSS LRG sample has led to the 
detection of acoustic oscillations in the spatial distribution of 
galaxies (Eisenstein et al. 2005), 23 it would be instructive to 
compare the "strength" of this survey to the planned blank sky 
SZ cluster surveys like PLANCK and SPT In Fig. 12 we show 
the number density of clusters as a function of redshift for the 
above described two types of SZ survey. The upper group of 
lines corresponds to the SPT-like deep survey while the lower 
curves are for a wide and shallow survey like PLANCK. For 
each of the surveys we have varied the lower mass cutoff M\ ovl 
so as to obtain in total 15,000, 25,000 and 35,000 clusters. 
These three cases are shown with solid lines. Dashed lines dis- 
play the pure flux-limited surveys without any lower cutoff in 
the mass imposed. Using these results and also taking into ac- 

23 The detection has also been claimed using the 2dF redshift survey 
(Cole et al. 2005). 
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Fig. 12. Comoving number density of clusters as a function of redshift 
for the PLANCK- (lower set of curves) and SPT-like (upper set of 
curves) SZ surveys. Solid lines in each set correspond to the cases with 
15, 000,25, 000 and 35, 000 clusters. The vertical dotted line shows the 
applied upper limiting redshift for the PLANCK. Dashed lines repre- 
sent purely flux-limited cases i.e. without any lower mass cuttoff im- 
posed. 



count the proper biasing factors as given by the square root 
of Eq. (18) divided by P(k, z = 0) 24 we can readily obtain 
AP/P as given by Eq. (6) and (7). The results of this calcula- 
tion are given in Fig. 13. Here the solid lines correspond to the 
PLANCK-like and dashed ones to the SPT-like surveys. Each 
set of lines corresponds to the detected cluster numbers (start- 
ing from above): 15, 000, 25, 000 and 35, 000. The dash-dotted 
curve, showing the results for the SDSS LRG sample, is found 
using Eq. (6) and the data for the effective volume given in 
Fig. 1 of Eisenstein et al. (2005) (again Ak = 0.005 /zMpc" 1 
was assumed). We can see that on large scales future SZ sur- 
veys have enough strength to improve the results obtained us- 
ing the SDSS LRG sample. However, the SDSS LRG sam- 
ple is going to double in size within a few years as the sur- 
vey is completed. The achievable AP/P for this final sam- 
ple is shown in Fig. 13 as a dotted line. Moreover, it seems 
that acoustic oscillations are able to survive at the quasilinear 
scales (k ~ 0.1 . . .0.3 /zMpc~') (see Fig. 4) which significantly 
increases the amount of information available for the galaxy 
redshift surveys. This is also confirmed by the recent N-body 
simulations by Springel et al. (2005) and Seo & Eisenstein 
(2005). In order to fully exploit this information one needs a 
complete theoretical understanding of how nonlinear effects, 
redshift space distortions and nonlinear biasing influence these 
features. So far there have been only a few works studying these 
important issues (e.g. Meiksin et al. 1999; Springel et al. 2005; 
Seo & Eisenstein 2005; White 2005) and we do not have a full 
theoretical description of them available yet. For this reason we 
have not attempted to incorporate the SDSS LRG sample into 
our Fisher matrix parameter estimation process. 

24 The light-cone bias parameters calculated this way are 5.0,4.7,4.3 
for the PLANCK and 4.4,4. 1,3.9 for the SPT with the number of clus- 
ters 15, 000,25, 000 and 35, 0000, respectively. 



12 



G. Hiitsi: Clustering of SZ clusters 



PLANCK 
SPT 
SDSS LRG 
SDSS2 LRG 




0.06 0.08 
k [h Mpc" 1 ] 



Fig. 13. Comparison of the "strength" of the future SZ cluster surveys 
with respect to the SDSS LRG. Lines shown for the PLANCK- and 
SPT-like surveys correspond to the detected cluster numbers (starting 
from above): 15, 000, 25, 000 and 35, 000. 



3.4. Some remarks on SZ vs. optical cluster selection 

SZ cluster selection might be superior to the simple optical 
one since SZ brightness does not suffer from ordinary cosmo- 
logical dimming. This allows one to obtain an approximately 
mass-limited sample of clusters that is spatially much more 
uniform than optically selected samples. For example, in Fig. 
12 the number density of clusters for both PLANCK (up to 
z ~ 0.6) and the SPT type of surveys drops approximately as 
oc z~ 2 5 . In contrast, the number density of the optically selected 
SDSS LRG sample (which contains mostly galaxies that popu- 
late dense cluster environments) drops as oc beyond z ~ 0.3. 
So in general with SZ-selected clusters one is able to probe 
larger volumes. Also as a larger part of the sample is at higher 
redshifts the clustering signal is stronger due to the increase of 
the bias factor with increasing distance. (Of course the ques- 
tion of whether these far away parts of the sample have high 
enough number density in order to be useful for the clustering 
study depends on the specific parameters of the experiment.) 
The other weaknesses of the optical selection are projection 
effects and confusion with the background objects. As an ex- 
ample, comparison of the X-ray and optically selected cluster 
catalogues often yields a rather poor match (e.g. Donahue et al. 
2002), which probably signals that many of the "optically con- 
structed" systems are actually false detections. 

For the SZ surveys the spectroscopic follow-up is a crucial 
issue, but several other studies (such as cluster number counts) 
also require determination of redshifts. As these investigations 
will be performed anyway, the clustering study can be seen as 
coming esentially "for free". 



4. Constraints on Dark Energy 

4.1. 2D power spectrum on a light-cone 

Having obtained a parametrized (and well calibrated) analyti- 
cal model for the light-cone power spectra of SZ-selected clus- 



ters of galaxies, we can estimate the accuracy with which it is 
possible to recover cosmological parameters. Since the obser- 
vations are done in cosmological (as opposed to the comoving) 
redshift space there are two additional effects one has to take 
into account: 

1. The increase of power along the line of sight due to the 
large scale coherent inflows towards massive accretion cen- 
ters. This effect is accounted for using results from linear 
theory. 

2. Cosmological distortion due to the fact that one can directly 
observe only redshifts, and in order to find the correspond- 
ing comoving distances, one has to assume some cosmo- 
logical model. Choosing an incorrect model will lead to 
distortions along and perpendicular to the line of sight 25 . 

Both of these effects will, in general, lead to an anisotropic 
power spectrum. Thus instead of a one dimensional 
(isotropized) power spectrum one has to consider here a two 
dimensional power spectrum with components along and per- 
pendicular to the line of sight. For simplicity, in the following 
we use a flat sky approximation. The 2D power spectra are cal- 
culated following the description given in Magira et al. (2000) 
with a slight modification to allow for various weight functions 
26 . Magira et al. (2000) took the fiducial cosmology given by 
the Milne's empty universe model. Thus the comoving distance 
intervals along and perpendicular to the line of sight are: 



Ax" 



ef 



Az, 



Ax? = — zA6, 



(25) 



(26) 



where Az is redshift interval and AO angular separation between 
two objects. For the general FRW universe the corresponding 
intervals read as: 

c 



Ax\\ = 



■Az, 



H(z) 
Ax ± = d M (z)A6, 



(27) 
(28) 



where c/m(z) = (1 +z)c/a(z) is the comoving transverse distance. 
Now, denning the shift parameters: 

Axii 

j(z), (29) 



qi(z) = 



A4 et 



(30) 



we can write down the final expression for the 2D power spec- 
trum on a light-cone: 

P LC ' 2D (k h k ± ;>M l0W ) = 



-max / , , w 

Jdzf W\z) l+/3(>M low ,z)(^J 



P c (k';> M low ,z) 



j dzf^W 2 (z)c ± (z) 2 C||(z) 



(31) 



25 The Alcock-Paczynski test (Alcock & Paczynski 1979) is based 
on these cosmological distortions. 

26 In the calculations we use the FKP weight function. 
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Fig. 14. Upper left panel: 2D light-cone power spectrum for clusters with mass above 1.0 ■ 10 14 /r'M o in logarithmic units. Lower left panel: 
as above, but the redshift distortion "switched off". Upper right panel: as upper left panel, only smooth component of the spectrum divided out 
and results shown using linear scale. Lower right panel: the same procedure applied to the lower left panel. 



where 

11 C||(Z)' ± < (-) V ii 

and 

P(> M low ,z) = -- 



1 



dlnD + (z) 



b eS (> M low ,z) dln(l+z) 



(32) 



(33) 



The factor c ± (z) 2 cy(z) in the denominator of Eq. (31) is the 
Jacobian determinant taking into account the change in a vol- 



ume element. It is missing in the numerator due to the cancella- 
tion by the similar but inverse term arising from the transforma- 
tion of the £-space volume element. The term in square brack- 
ets models the amplification due to the coherent inflows and the 
last term in the numerator, P c (k'\ > M\ ovl ,z), is given earlier by 
Eq. (20). In the case of SZ flux-selected clusters, M\ ov/ at each 
redshift for a given lower flux limit is found using Eq.(22) and 
(24). 

In Fig. 14 we present some examples of 2D power spec- 
tra calculated in the manner described above. The top left- 
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Fig. 15. Error ellipses for the PLANCK-like SZ survey including CMB priors as described in the text. Solid/dashed lines show results for 3/1 
redshift bin(s). 



hand panel shows the 2D light-cone power spectrum of clus- 
ters with a mass above 1.0 • 10 14 /z _1 M G (in logarithmic units) 
while the lower left-hand panel contains the same spectrum 
but with the redshift space distortion "switched off". The con- 
tours starting from the upper right corner correspond to the 
values 10 4 ' 8 and 10 4 7 for the upper and lower panel, respec- 
tively and the step size was taken 10° 1 h~ 3 Mpc . It is clearly 
seen how linear redshift-space distortion boosts power along 
the line of sight. The cosmological distortion (in the currently 
selected reference model) on the other hand works in the op- 
posite way. Since the chosen Milne model is strongly different 
from the ACDM "concordance" model the cosmological dis- 
tortion is easily visible. We use Milne's model only for illus- 
trative purposes but in the following parameter estimation part 
we change the reference model to the WMAP "concordance" 
cosmology. On the right-hand panels of the figure we have re- 
moved the smooth component of the power spectrum revealing 
the series of acoustic rings. Here the continuous contours cor- 
respond to P LC ' 2D /P^^ h values of 0.95 and 1.05 while the 
dashed lines are for the values 0.97 and 1.03. It is important to 
note that the picture on both panels looks practically the same. 
This is due to the fact that cosmological and redshift-space dis- 
tortions work in a different way: cosmological transformation 
stretches or compresses the power spectra on the plane of the 



figure whereas redshift distortion moves the spectra in a verti- 
cal direction. If we had smooth power spectra without any par- 
ticular features then it would be extremely hard to disentangle 
these two types of distortions. Having the power spectra with 
acoustic features it is easy to isolate cosmological distortion by 
dividing out a smooth component. The ability to disentangle 
cosmological and redshift-space distortions is extremely im- 
portant to extract the bias parameter from the survey in a self- 
consistent way. In the parameter estimation part of this section 
we see how much better one does with the model having acous- 
tic oscillations as compared to the one without. 



The total power spectrum measured over a broad z-interval 
is a weighted sum of differently distorted power spectra and 
so there will be some loss of acoustic features. The loss is 
stronger along the line of sight, as can be seen from Fig. 14. 
This is due to the currently chosen reference model where 
\—jir\ > l^jpl- Again, this effect is strongly pronounced be- 
cause Milne's model differs strongly from the WMAP "concor- 
dance" cosmology for which our calculation was done. 
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Fig. 16. Analog of Fig 15 for the SPT-like survey. Solid/dashed lines show results for 4/ 1 redshift bin(s). 



4.2. Parameter estimation 

In this subsection we apply a Fisher matrix forecasting tech- 
niques to study how well one can determine cosmological pa- 
rameters. Since the pioneering investigations in the field of 
CMB anisotropics (Jungman et al. 1996) and galaxy redshift 
surveys (Tegmark 1997) these methods have gained great popu- 
larity. For a full description of the method with applications see 
Tegmark et al. (1997). Following Tegmark (1997) the Fisher 
matrix in the case of a 2D power spectrum can be written as 
(see also Hu & Haiman (2003)): 



^ = ZZ 



dlnP LC - 2D (fc|| n ,fc ±m ) 



<LC,2D/ 



||n> ^±m) 



where 



(34) 



(35) 



and we take Ak ± = Afey = Afe. p LC > 2D is given in Eq. (31) and 
V e ff in Eq. (7). The wavevector components fey; and k±j form 
a rectangular grid with a step size Afe. We allow them to span 
the range 0.005 . . .0.1 /zMpc" 1 and take Afe = 0.005 /zMpc" 1 



i.e. we have a 20x20 grid. It is not justified to go to higher 
wavenumbers than 0.1 /zMpc -1 since there already the simple 
linear scaling of the cluster power spectrum seems to break 
down (see Fig. 6). The parameters ; are taken to be £2 m , cr 8 , 
Q h , h, wo plus the bias parameters b n in each redshift bin. Thus 
in the case of four redshift bins (which is the maximum num- 
ber considered in our analysis) we have a total of nine free pa- 
rameters. We also assume a fiducial cosmology given by the 
best-fit WMAP "concordance" model (Spergel et al. 2003) plus 
the dark energy in the form of the cosmological constant i.e. 
w = -1. 

Here we perform calculations for two types of SZ survey. 
The first one that should serve as a prototype for PLANCK is 
a shallow survey that covers the full sky with a sensitivity limit 
of Fi ow = 30 mJy at 353 GHz 11 . We also apply a lower mass 
cutoff of Mi ow = 2.75 • l0 u h~ l M Q resulting in total ~ 25,000 
clusters up to the applied limiting redshift z - 0.6 (see Fig. 
10 lower panel). The calculations are done for one redshift bin 
spanning z = . . . 0.58 and also for three bins: z\ = . . . 0.33, 
z 2 = 0.33 . . .0.45, z 3 = 0.45 . . .0.58. The achieved AP/P val- 
ues at fe = 0.05 /zMpc -1 are 4.3% for the one bin case and 
7.5% for each of the bins in the three bin case. The other 
deep and rather "narrow" survey should mimic the performance 



http://astro.estec.esa.nl/Planck 
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Fig. 17. Upper panel: Error ellipses under various assumptions for the 
PLANCK-like survey. Here we have shown the following cases: 1. 
Clustering signal of the model without baryonic oscillations + CMB 
priors, 2. #1 + prior on bias, 3. Clustering signal of the "wiggly" 
model, 4. #3 + CMB priors, 5. #4 + prior on bias. Lower panel: as 
above, except for the SPT-like survey and with a changed order of the 
cases #2 and #3. Lower panel inset: constraints obtainable with: 1. 
SPT, 2. PLANCK, 3. SPT + PLANCK, 4. The survey with SPT char- 
acteristics covering the full sky. All of the cases here have CMB and 
bias priors added. 



of SPT. Here we have taken F\ ow = 5 mJy at 150 GHz and 
Mi ow = 1.5 • IO'VMq (see Ruhl et al. (2004); Majumdar & 
Mohr (2004)). Moreover, we assumed that the survey is capa- 
ble of covering one octant of the sky and here we have not ap- 
plied any high redshift cutoff. As for the PLANCK-like survey 
two different cases are considered here. In the first case while 
combining all the data to one common power spectrum we are 
able to achieve AP/P ^ 5.1%, whereas in the second case 
with the four redshift bins: z\ = . . . 0.45, zi = 0.45 . . . 0.63, 
z 3 = 0.63 . . . 0.83, za = 0.83 ... we find a AP/P =* 9.9% in each 
redshift bin. The total number of objects detected with this type 
of survey would be ~ 27, 000. 

The results in Fig. 15 and 16 present constraints for the five 
cosmological parameters: Q m , cr 8 , £2 b , h and wq. Here we have 
added priors to O m /i 2 and Q.\,h 2 from the CMB angular power 



spectrum studies, and also a prior to the bias parameters. In 
order to be able to compare our results to the ones given in 
Hu & Haiman (2003) the fractional errors of 0.01 for O m /; 2 
and Ob^ 2 (this should be achievable with the PLANCK mis- 
sion (Hu 2002)) were similarly assumed. These constraints are 
easily "rotated" to the frame used here since under coordinate 
transformations the Fisher matrix transforms as a second rank 
tensor. Moreover, we have restricted our calculations to the flat 
models only, and for the bias have assumed that one is able to 
describe it with a relative accuracy of 15%. Finally, the joint 
Fisher matrix is the sum of all the Fisher matrices transformed 
to a common frame. All the calculations done here assume an 
underlying model with baryonic features in the matter power 
spectrum. Fig. 15 presents results for the previously described 
shallow survey with full sky coverage while the error ellipses 
in Fig. 16 apply to the deep and narrow survey. In both fig- 
ures dashed lines correspond to the single and solid lines to the 
multiple bin case. Due to the fact that with a single bin one is 
able to measure the shape of the power spectrum with a higher 
precision than in the case of multiple bins, we see from the 
above figures that in general stronger constraints on O m are ob- 
tained. On the other hand constraints on wo are much stronger 
in the multiple bin case due to the increased knowledge about 
the redshift derivatives. We also performed calculations taking 
two redshift bins and the results were already rather close to 
the solid curves in the figures above which can be interpreted 
as an indication that any further redshift slicing would not im- 
prove constraints on wq. Also one should not increase the num- 
ber of redshift bins much above the maximally used values 
of three and four since then the wavevector bins with width 
Ak = 0.005 hMpc 1 would become highly correlated and the 
above Fisher matrix calculation would not be meaningful. 

Probably the most interesting constraints in Fig. 15 and 16 
are the ones for Q m (or Qqe since we have assumed flat models) 
and wo. Error ellipses in the Q m -wo plane are also given in Fig. 
17 for various different assumptions. The top panel here cor- 
responds to the PLANCK-type and the lower one to the SPT- 
type of survey. The order of ellipses in the top panel starting 
from the bottommost one is as follows: (1) clustering signal of 
the model without baryonic oscillations + CMB priors, (2) (1) 
+ prior on bias, (3) clustering signal of the "wiggly" model, 
(4) (3) + CMB priors, (5) (4) + prior on bias. The only dif- 
ference in the lower panel is the reversed order of (2) and (3). 
The inset in the lower plot displays the constraints obtainable 
(again starting from the bottommost ellipse) with (1) SPT, (2) 
PLANCK, (3) SPT + PLANCK, (4) the survey with SPT char- 
acteristics covering the full sky. For all of the cases shown in 
the inset we have included CMB and bias priors. The constraint 
ellipses for the "clustering only" case assuming a "smooth" 
model would fill almost all the plot area and for the sake of clar- 
ity we have not displayed them here. It is evident from Fig. 17 
that the model with baryonic oscillations is performing much 
better compared to its smoothed counterpart. Adding prior in- 
formation (in contrast to the "smooth" case) does not result 
here in a strong improvement i.e. the clustering signal alone 
already has a significant constraining power. 

Table 1 lists the principal components of the clustering 
analysis only i.e. no CMB and bias priors included. Also we 
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Table 1. Principal components as given in Eq. (36) for the "clustering only" case for PLANCK-like, SPT-like and for the combined survey. 





i 


X 


e \ 


e 2 


e\ 


e\ 


e s 




1 


1.6 


10 3 


0.74 


-0.0064 


-0.24 


0.62 


0.090 


PLANCK 


2 


1.7 


10 2 


-0.48 


-0.015 


-0.43 


0.49 


-0.58 


(3 z-bins) 


3 


32. 




0.38 


-0.53 


-0.23 


-0.46 


-0.56 




4 


18. 




0.027 


0.41 


-0.80 


-0.37 


0.26 




5 


5.9 




0.27 


0.74 


0.28 


-0.13 


-0.53 




1 


9.5 


10 2 


0.71 


-0.094 


-0.23 


0.63 


0.21 


SPT 


2 


85. 




0.15 


-0.088 


0.29 


-0.36 


0.87 


(4 z-bins) 


3 


54. 




-0.012 


-0.99 


-0.032 


-0.10 


-0.13 




4 


17. 




0.69 


0.094 


0.24 


-0.54 


-0.41 




5 


4.6 




0.049 


0.057 


-0.90 


-0.42 


0.12 




1 


2.6 


10* 


0.73 


-0.038 


-0.24 


0.63 


0.13 


PLANCK+ 


2 


2.5 


10 2 


-0.38 


0.044 


-0.39 


0.45 


-0.70 


SPT 


3 


78. 




0.18 


-0.92 


-0.027 


-0.22 


-0.28 




4 


52. 




0.54 


0.39 


0.087 


-0.44 


-0.60 




5 


24. 




0.021 


0.057 


-0.89 


-0.40 


0.22 



have marginalized over bin bias parameters. This serves as a 
compact way of summarizing our results. The principal com- 
ponents are given in the form: 



SPT+PLANCK case, however, in the case of Q m an approxi- 
mate agreement is achieved with our SPT full sky example. 



n 

j=i 



©fid 



= i + 



i 



for all 



1 ...5. 



(36) 



Here ej is the j-th component of the z'-th eigenvector and 



A 1 is the ;-th eigenvalue of the matrix F\\ 



@fid fld F 



(no summation over indices). The parameter vectors 
(Q m ,cr 8 ,Q b ,/z,w ) and fid = (0.27,0.84,0.044,0.71,-1.0). 

4.2.1 . Comparison to previous work 

Probably the two closest works to ours are Majumdar & Mohr 
(2004) and Hu & Haiman (2003). In Majumdar & Mohr (2004) 
the authors discuss constraints obtainable by combining the 
cluster power spectrum with independent information from 
cluster number counts. Unfortunately they do not present re- 
sults separately for the clustering signal only. Moreover, they 
use an isotropized power spectrum which leads to a signifi- 
cant loss of information, especially when the spectra contain 
baryonic features. In Hu & Haiman (2003) the authors use the 
full 2D power spectrum of galaxy clusters although they do not 
take into account light-cone effects. Since the light-cone power 
spectrum is a blend of differently deformed power spectra (if 
we are away from the reference model point) some loss of bary- 
onic features will result as seen from Fig. 14. Also in their anal- 
ysis they used fc-modes up to 0.15 as opposed to our adopted 
value of 0.1 /zMpc -1 . As seen from Fig. 6 at k = 0.15 /zMpc -1 
the cluster power spectrum already differs quite significantly 
from the simple linear one, particularly for the more massive 
systems. To simplify the comparison our results with Hu & 
Haiman (2003), we used identical CMB priors. An additional 
difference is that we do not allow the spectral index of the 
power spectrum to vary but keep it fixed to n = 1. In total 
their results for the SPT-type of survey are significantly more 
optimistic, e.g. the constraints on wq agree roughly with our 



5. Conclusions 

In this paper we studied the clustering of SZ-selected galaxy 
clusters on a past light-cone with particular emphasis on con- 
straining the properties of DE. We implemented an extended 
Press-Schechter type of analytical model as described e.g. in 
Sheth & Tormen (1999). The description of the calculation 
of the light-cone power spectra (e.g. Yamamoto et al. 1999) 
was modified slightly to incorporate other than simple num- 
ber weighting schemes. The analytical model was extensively 
calibrated using the outputs from the VIRGO Consortium's 
Hubble Volume simulations. With a little bit of fine tuning we 
were able to match analytical light-cone power spectra with 
ones extracted from the simulations to an accuracy better than 
20%. The SZ scaling relations were calibrated so as to get a 
good match to the number count results from the state-of-the- 
art hydrodynamical simulations of White et al. (2002). Having 
a well calibrated analytical model we investigated how accu- 
rately future SZ surveys like PLANCK and SPT could deter- 
mine the cluster power spectrum and whether they would be 
able to detect traces of baryonic oscillations. Also we made 
use of VIRGO simulation outputs to build cluster catalogs for 
various survey depths and sensitivity limits. We showed that 
the aforementioned blank sky SZ surveys will be able to im- 
prove the detection of acoustic features based on the SDSS 
LRG sample. To obtain a high-fidelity detection of the bary- 
onic oscillations, one has to reach a relative accuracy ~ 5% if 
the wavenumber bin width Ak = 0.005 ZzMpc -1 is used. This 
seems hard to achieve with surveys like SPT still having rela- 
tively narrow sky coverage. On the other hand, for surveys with 
a wide sky coverage the prospects seem to be better e.g. with 
~ 25, 000 most massive clusters up to redshift ~ 0.6 one should 
be able to reach a relative accuracy ~ 4.5% at k — 0.05 h Mpc~' 
i.e. roughly the scale where one expects to see the first major 
acoustic feature. 
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In the last part of the paper we carried out a Fisher matrix 
forecasting analysis for cosmological parameters, which add up 
to nine in the case of fourredshift bins: £2 m , <x 8 , Q.^,, h, wo plus a 
free bias parameter for each bin. We included prior information 
for £l m h 2 and Q. h h 2 from CMB angular power spectrum studies 
and also constrained the possible values for the bias parame- 
ters. A prior on bias parameters only has a significant effect 
in the case of models with smooth power spectra i.e. models 
with acoustic oscillations have enough constraining power to 
give an estimate of bias parameters from the survey itself. The 
most interesting constraints are obtained for Q m and wq. Wide 
and rather shallow surveys like PLANCK in combination with 
a CMB prior on Q. m h 2 are able to provide strong constraints on 
Q. m or in the case of flat models equivalently on Q DE = 1 - £l m . 
The constraints on wq on the other hand are not as good as the 
ones obtained by deeper and narrower surveys with the charac- 
teristics of SPT due to the lack of higher redshift objects. We 
also give our results as the principal components of the Fisher 
matrix that should allow for an easy way of comparison and 
also for a fast way of incorporating these constraints into fur- 
ther parameter forecasting studies. 
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